Electron energy loss in carbon nanostructures 

A. Rivacoba 1 ' 2 and F. J. Garcia de Abajo 2 
1 Matenalen Fisika Saila, Kimika Fakultatea, UPV/EHU, 1072 P.K., 20080 Donostia, Spain 
2 Donostia International Physics Center (DIPC) and Centro Mixto CSIC-UPV/EHU, Aptdo. 1072, 20080 San Sebastian, 

Spain 
(February 1, 2008) 

The response of fullerenes and carbon nanotubes is investigated by representing each carbon atom 
by its atomic polarizability. The polarization of each carbon atom produces an induced dipole 
that is the result of the interaction with a given external field plus the mutual interaction among 
carbon atoms. The polarizability is obtained from the dielectric function of graphite after invoking 
the Clausius-Mossotti relation. This formalism is applied to the simulation of electron energy loss 
spectra both in fullerenes and in carbon nanotubes. The case of broad electron beams is considered 
and the loss probability is analyzed in detail as a function of the electron deflection angle within 
a fully quantum-mechanical description of the electrons. A general good agreement with available 
experiments is obtained in a wide range of probe energies between 1 keV and 60 keV. 
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I. INTRODUCTION 

Since the discovery of fullerenes, 1 the collective exci- 
tations of these molecules have received considerable at- 
tention. The first electron energy loss spectra (EELS) 
of solid Ceo where reported by Saito and coworkers. 2 
Later, EELS has been employed to study fullerenes sup- 
ported on surfaces, 3:4 single-walled fullerenes in the gas 
phase, 5 multishcll fullerenes, 6-9 and other carbon struc- 
tures of more exotic shapes. 10-15 Along with this ex- 
perimental effort, theoretical studies of the electronic 
structure and the collective excitations of these molecules 
have also been carried out. 16-19 In particular, Saito and 
coworkers 2 noted the strong similarity between loss spec- 
tra of fullerenes and graphite, and pointed out that 
this was due to the fact that the electronic structure 
of both of these materials consisted of a orbitals along 
the carbon bonds and tt orbitals in the normal direc- 
tion. This idea has been exploited by Lucas et al. e 
to describe the response of fullerenes in connection to 
EELS. They used a continuum model consisting of a clas- 
sical dielectric sphere that was made of a non-isotropic 
dielectric function equal to that of graphite. Henrard 
and Lambin 20 extended that work by using the discrete- 
dipole approximation 21,22 (DDA), where the fullerenes 
were considered as a system of coupled point dipoles, and 
their atomic polarizability tensor was obtained from the 
dielectric function of graphite via the Clausius-Mossotti 
relation. 

The continuum model 6 works well for multishell 
fullerenes, where the thickness of the target can be prop- 
erly assigned, but this is not the case of single-layer 
clusters, where the thickness parameter becomes critical 
while it is not well defined. However, the DDA provides 
a reasonable description in the latter case. 

So far, EELS electrons have been described as classical 
charged particles, 20 although such an approach presents 
some disadvantages, connected to the fact that the en- 
ergy loss spectra depend strongly on the beam position 



relative to the specific atomic positions of the carbon 
cluster, and in particular, the loss probability diverges 
when the beam intersects any of the atomic positions, a 
fact that derives from the unphysical singularity at the 
origin of the dipole potential (~ 1/r 2 ) used in this ap- 
proach. This pathological behavior at small distances is 
connected to unphysically large values of the momentum 
transfer, which are included implicitly in the classical en- 
ergy loss theory 23,24 The classical theory is then only jus- 
tified for relatively distant collisions, where the contribu- 
tion of large momentum transfer is negligible. Another 
drawback of a classical description of the probe derives 
from its inability to account for the momentum transfer 
dependence of the loss probability, a magnitude which 
is accessible to EELS theory 25 and experiments 7,9 ' 15 via 
the selection of electron scattering angles: momentum- 
resolved EELS has actually proved to be a useful tool 
for obtaining optical and electronic properties in many 
cases. ' ' One more issue related to momentum transfer 
(in this case with the component parallel to the electron 
velocity) is the so-called recoil of the electron, which can 
play a noticeable role at low electron energies like those 
employed in some experiments with fullerenes in order to 
avoid damage of the specimen. 23 

In this work, we make use of the self-energy formal- 
ism to calculate the energy loss probability within a fully 
quantum-mechanical description of the electron, while 
the response of the carbon cluster is expressed in terms 
of the screened interaction as calculated in the DDA via 
the atomic polarizability of the carbon atoms. The en- 
ergy loss probability is obtained as a function of both 
the lost energy and the electron deflection angle (i.e., the 
momentum transfer component normal to the electron 
beam). This solves the problems that come about in a 
classical description of the electron, as discussed above: 
divergences in close encounters are avoided in a natu- 
ral way and the contribution of electron recoil is fully 
accounted for. For the electron energies under considera- 
tion, the initial and final electron states can be described 
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by plane waves. This is a very realistic approach to EELS 
experiments when a broad beam is employed. Further- 
more, when narrow beams of ~ 0.5 nm are used, one 
can interpret our extended-beam results as spectra aver- 
aged over the impact parameter, 23 which raises the ques- 
tion of the suitability of the dipole potential ~ 1/r 2 near 
the atoms; nevertheless close-encounter processes corre- 
spond to large momentum transfers, and therefore, they 
contribute little when calculating the spectra of electrons 
collected with small scattering angle, as it is usually done 
in the experiments. Notice that for classical electrons 
the inelastic scattering occurs up to an impact parame- 
ter b < v/u>, where v is the electron velocity and ui is the 
energy loss, so that the broad beam approach is likely 
to overestimate the weight of low-energy excitations for 
target sizes comparable to the lateral coherence length of 
the electron beam. 

The screened interaction calculated in Sec. II is actu- 
ally the Green function of Poisson's equation in the pres- 
ence of the sample, and thus, it provides the response 
of the fullerene to any arbitrary external field. When 
we apply this theory to simulate EELS in single layered 
carbon structures (Sec. V), very good agreement with 
experiments is obtained, indicating that the screened in- 
teraction could equally be applied to other problems in- 
volving the response of fullerenes with inclusion of atomic 
details (e.g., in STM, light absorption, 26 van der Waals 
forces, 27 and image states 28 ) . A similar approach could 
be also applied to other materials such as boron-nitride 
structures. 

Atomic units (a.u., i.e., U — m = e = 1) will be used 
throughout this work, unless otherwise specified. 



Notice that matrix F is fully symmetric both in coordi- 
nate indices ij and in atomic ones ab. From Eq. (1), the 
dipolar moments of the system are found to be 

p?M=2[M- 1 ]?/a5 fc £7 fc (r i , W ) ) (3) 
b,jk 

where M _1 is the inverse of the w-dependent matrix 
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Eq. (3) permits calculating the induced-dipolc mo- 
ments for an arbitrary external field, from which one can 
in turn obtain the induced field. In particular, this per- 
mits to derive the so-called screened interaction of the 
system W(t,t' ,u>), which is the solution of the Poisson 
equation for a point charge at r', that is, 



V 2 W(r,r',u) = -47r<5(r-r'). 



(5) 



In this case, the field E(r) created by the external 
probe is just the Coulomb field 1/ | r — r' |, and there- 
fore, the induced part of the screened interaction can be 
straightforwardly written as 
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II. THE SCREENED INTERACTION OF A 
FINITE ATOMIC CLUSTER 

We shall consider a system composed by N atoms lo- 
cated at the positions r a , with a = 1, . . . N. Each atom 
can be polarized so that it behaves like an induced dipole. 
The polarizability tensor of atom a will be denoted a a (co) 
in frequency space lu. When an external electric field 
E(a>) is applied to the system, the induced dipole mo- 
ment p a (uj) of atom a is the response to the total field, 
that is, the external field plus the field induced by the 
rest of the dipoles. One finds, 

jk 

where the indices i, j, and k refer to Cartesian coordi- 
nates. The second term in the right hand side of Eq. (1) 
involves the dipole-dipole interaction, which can be writ- 
ten in terms of the relative atomic positions r ab = r a — r b 
as 



F a t = — - — [ 
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a^b, 
a = b. 



Like M above, this matrix is symmetric in both sets 
of indices, so that W lnd (r, r', to) is also symmetric with 
respect to the spatial variables r and r'. 

From a practical point of view, it is important to note 
that the dependence on frequency and spatial variables 
is separated in this formalism: the dependence on lu is 
fully contained in matrix C, whereas the dependence on 
both r and r' comes from the last two factors in Eq. 
(6). Therefore, although the actual dipoles have to be 
calculated numerically to obtain matrix C, the spatial- 
dependent part can be dealt with analytically in many 
cases. For instance, it will be useful to rewrite (6) in 
momentum space as 
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The response of the system to any external probe can 
be now expressed in terms of Eq. (6). As an example, the 
energy loss experienced by a fast classical electron is ob- 
tained as a double integral of the screened interaction, 2 
which can be solved analytically 29 to recover expressions 
(18)-(20) of Ref. [ 20]. 
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III. ENERGY LOSS OF FAST ELECTRONS IN A 
SYSTEM OF COUPLED DIPOLES 

As an application of the formalism presented above, 
we now proceed to compute the energy loss experienced 
by a broad beam of fast electrons incoming along the z- 
axis into an atomic cluster, described as a set of coupled 
induced dipoles. We study the transition between ini- 
tial and final states of the electrons, described by plane 
waves * (r) = y-VV k ° r and */(r) = V^e^f, re- 
spectively, where ko and k/ are the initial and final mo- 
menta, and V is the normalization volume. The transi- 
tion probability can be written 32 ' 33 

P H = ^E / drdr'*?(r , )*/(r) 
/ J 

x Im{-W r i nd (r,r / ,w)}*S(r)*o(r / )%-wo/), (9) 

where v = ko is the velocity of the incoming probe, taken 
to be along the positive z axis, Uof = ^(k 2 , — k 2 ) is the 
transition energy, and Im{W md (r, r', u>)} is the imagi- 
nary part of the screened interaction. The sum over final 
states / can be transformed into an integral over the mo- 
mentum transfer Q = k — k/ . One finds 



iQ(r-r') 



x Im{TF nd (r, r', u)} 5(io - Qv + -Q 2 ) . (10) 
Now, inserting Eq. (8) into this expression, one can write 
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x J dQ 



^ e -^ 6 - ra )^-Qv + iQ 2 ). (11) 



The probability of collecting an electron after suffer- 
ing an energy loss u> and being scattered inside a circular 
aperture of half-angle 6 m can be expressed as the sum of 
the contributions of excitations characterized by w and q 
(the energy and the component of the momentum trans- 
fer normal to the beam direction, respectively, that is, 
Q = (l> Qz)) as follows: 



(12) 
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where Tjl is a matrix (symmetric in both ab and ij) that 
reduces to 
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dq 
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when a = b. Here, the z component of the momentum 
transfer, q z , and A are both functions of u> and q, that 
is, 



= v — \Jv 2 — 2uj 



A = (^ + aVl- 



2cu + q 2 



(14) 
(15) 



and q m ~ v9 m is the largest possible value of q. The 
matrix elements coupling different atoms can be written 
as 
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where J n is the Bcssel functions of first kind and order n, 
Pab = V(x a - x b ) 2 + {y a - y b ) 2 , and 6 ab = co S - 1 {(x b - 
x a )/pab}- 

The expression (14) for q z derives from the 5— function 
ensuring energy conservation in Eq. (11), where the term 
\Q 2 represents the recoil of the electron. For large elec- 
tron velocities, this term is negligible and q z can be ap- 
proximated by the classical value q z = ui/v. Actually, 
Ritchie 23 proved that even when all the scattered elec- 
trons were collected the effect of this term is small for 
fast (100 keV) electrons. 

The coupling between different atoms in Eq. (12) is 
introduced via the screened interaction, where the polar- 
ization of neighboring atoms is strongly correlated. This 
leads to the existence of polarization waves that play the 
same role as plasmons in extended media. The spatial 
extension of these waves in the direction normal to the 
beam can be analyzed by studying the dependence of 
the T^j integrals (see above) on the separation between 
dipoles in the direction normal to the electron beam, p ab , 
and for the sake of obtaining an estimate one can use the 
classical-approach approximations consisting in neglect- 
ing the recoil (i.e., q z = oj/v) and extending the integral 
in q from zero to infinity for a ^ b; one finds 
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where Kq is the modified Bessel function, and similar ex- 
pressions involving K m are found for the rest of the 
integrals. Taking into account the asymptotic behaviour 
of K m , Eq. (16) permits us to say that the dipoles are 
effectively coupled up to a distance of the order of the 
adiabatic length v/uj, similar to the lateral extension of 
the charge density associated to surface plasmons in ex- 
tended media. 30 

Equation (12) can be applied to any atomic cluster 
described by the atomic polarizabilities. In the case of a 
system of non-interacting dipoles [i.e, when is set to 
in Eq. (4)], assuming an isotropic polarizability, a" fe = 
a Sjk , and neglecting the recoil, the energy loss per atom 
probability reduces to 



(17) 



This expression reproduces the energy loss probability 
of a system of uncoupled dielectric spheres of radius R, 
where the polarizability is given by a — fq^-R 3 , in the 
limit of small momentum transfer, qR -C 1. 31 > 32 > 34 This 
means that the suitability of this approach is restricted 
to values q m < 1/R, where R has to be understood as 
the radius of a region containing the valence electrons. 
For a electrons one can estimate R as half the distance 
between C atoms in graphite (i.e., 1.42 A), which limits 
the validity of the present theory to q m < 1 a.u. No- 
tice that this limit is related only to the suitability of 
the dipolar potential near the origin: as one reaches the 
limit qR ~ 1 , higher multipolar excitations become more 
relevant than the dipolar one, and one should rely on a 
multipole approach. 35 This fact, pointed out by Keller 
and Coplan, 5 posses another limit to the applicability of 
the DDA theory. 

Although it has been proved 24 that classical energy loss 
expressions can be recovered from quantum-mechanical 
ones by neglecting recoil and taking the q m — > oo limit, 
such an expression is not defined in the present context, 
since diagonal terms of T diverge logarithmically at large 
values of q, which is another manifestation of the unphys- 
ical behaviour of the dipolar potential near the target- 
atom nuclei. 



IV. THE SCREENED INTERACTION OF 
INFINITE SYSTEMS: CARBON NANOTUBES 
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where (j>f b (q) is the Fourier transform of F^ b , which re- 
duces to 
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with i,j = 2,3 (j = 1 refers to the x direction) and 
Pab = \/(y a - y b ) 2 + {z a - z b ) 2 . Notice that <j>$(q) de- 
pends only on the components of the relative position 
vector that are normal to the translational axis of sym- 
metry, and therefore, it is enough to compute it for the 
atoms of a given unit cell of the tube. 

For atoms located in equivalent positions of the tube 
(i.e., p a b = 0), only the diagonal elements are involved, 
but obviously, the above expressions diverge as 
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However, this divergence can be avoided on physical 
grounds by excluding the self- interaction (i.e., an atom 
interaction with itself) from the summation over target 
atoms [see Eq. (2)]. In a similar way, one can write 



[m- 1 ^ = r 

Jo 



dqe lq{xa - xb) v? b (q) 



(21) 



where g = 2n/d is the reciprocal lattice constant and 
v?j{q) satisfies the secular equation 

- hi S ac = £{ % S ab < fe Rt b (q)}v b f(q) . (22) 
9 b,kj 

Here, the indices a, 6, and c run only over the first unit 
cell and R^ b (q) is defined as 



> ing(x a -x ) lab 



(23) 



The theory presented in the preceding sections can be 
straightforwardly applied to finite atomic clusters, such 
as fullerencs, but it is not suitable to study infinite sys- 
tems. In what follows, we will recast this theory to de- 
scribe a carbon nanotube consisting of the infinite repe- 
tition of a unit cell of length d along the direction of the 
x-axis. In order to take advantage of the translational 
invariance of this system, we write the coupling integrals 
F*» as 



For a = b (equivalent atoms), the n — term must be ex- 
cluded from Eq. (23) in order to avoid the self-interaction. 
This removes the noted divergence from Eq. (20) and al- 
lows using finite numbers in the secular equation (22). 
Finally, one finds 



RTM - - -UK® = ± £ 



n>0 



(24) 
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which is a well defined sum. The fact that R$(q) is a 
periodic function of period g has been exploited to re- 
duce the integration interval in Eq. (18) to the first Bril- 
louin zone in Eq. (21). Once the matrix v(q) has been 
obtained, M" 1 can be computed via Eq. (21) and the 
problem is solved. 

V. RESULTS AND DISCUSSION 

In order to apply this theory to fullerenes, one must 
obtain first the polarizability of the C atoms in the 
molecule. This is done via the Clausius-Mosotti rela- 
tion by using the bulk dielectric function of graphite, 36 
assuming that the polarizability of C atoms is essentially 
the same in graphite and in fullerenes owing to their sim- 
ilar local atomic environment. In Fig. 1, we plot both 
the real and the imaginary parts of the two indepen- 
dent components of the polarizability tensor (i.e., normal 
and parallel to the graphite planes). These plots corre- 
spond to C atoms in a graphite-like structure and have 
no physical meaning unless the coupling with neighbor- 
ing atoms is considered. It is interesting to note that 
both components present resonances around 6 eV, which 
lead to the excitation usually referred to as tt, although 
the coupling with the a electrons in the plane is also im- 
portant, as shown by the actual position of this peak in 
the EELS data discussed below. On the other hand, the 
broad peak of the imaginary part of the parallel compo- 
nent around 27 eV reproduces correctly the loss spectrum 
of graphite. 2 This correspondence between some polariz- 
ability features and EELS experiments reflects the fact 
that a electrons are almost uncoupled at energies above 
10 eV in graphite: the normal component of the polar- 
ization is very small and the coupling between a elec- 
trons within a carbon plane vanishes due to the symme- 
try of the lattice (the coupling between atoms in different 
planes is also very weak due to the large separation be- 
tween planes). Nevertheless, the surface curvature in car- 
bon nanostructures breaks that symmetry, and therefore, 
the spectra become sensitive to the geometrical shape of 
the sample over distances larger than the interatomic sep- 
aration. 

Let us first study the energy loss of an electron beam 
in a (10, 10) carbon nanotube to illustrate a case of an 
infinite atomic structure. 19,37 The basic cell of this tube 
consists of a ring of diameter 1.37 nm, formed by 40 C 
atoms. This ring is periodically repeated with a period of 
d = 0.247 nm along the tube axis. For each C atom, the 
polarizability tensor has been built so that its normal 
axis is perpendicular to the plane formed by the three 
closest atoms. We have considered a typical experiment 
in transmission electron microscope (TEM), where a 100- 
keV electron beam is impinging normal to the cylinder 
axis. The screened interaction has been calculated in 
momentum space as shown in Sec. IV, and we work out 
the energy loss probability per C atom by taking into ac- 



count the contribution of neighboring cells to the sum of 
expression (12) [ten cells at each side of the central cell 
are enough to get good convergence above 3 eV (see Fig. 
2(b))]. Fig. 2(a) shows the energy loss probability for 
a typical value of the collector aperture of 9 C = 1 mrad. 
The spectrum exhibits three features around 6.5, 11, and 
between 17 — 20 eV, which according to the analysis of 
Rcf. [ 14] can be identified as 7r-surface plasmons and 
both low- and high-energy a + tt surface plasmons, re- 
spectively. Below the tt feature, there is a shoulder in 
the loss spectrum that has been experimentally observed 
in both fullerenes 3-5 ' 10 ' 11 and nanotubes. 14 ' 15 Compar- 
ison with Fig. 1 permits to postulate that these peaks 
derive from the coupling between the resonances of the 
tt and a electrons. 

In Fig. 2(a), we have also plotted the contributions 
to the sum of Eq. (12) of the first neighboring cells with 
respect to a given central cell. Although the coupling be- 
tween different cells involves atoms at different distances 
from each individual atom, these plots provide valuable 
information on the distance dependence of the electro- 
dynamic interaction between atoms. First, one can see 
that the central cell and the first neighboring cells are 
give contributions of the same order of magnitude. The 
strong cancellation between contributions from the cen- 
tral cell and the first neighbors at u> = 5 eV (i.e., the 
position of the sharp resonance in the normal component 
of the polarizability) proves the collective nature of tt- 
plasmon resonances. This plot shows that the coupling 
between nearest cells gives rise almost entirely the to- 
tal energy loss probability between 7 eV and 17 eV, but 
more distant cells are needed outside this energy range. 
This point is relevant to the study of the a + tt plasmon, 
which is around 24 eV in fullerenes, whereas it shows up 
around 16 eV in nanotubes. As discussed above, this en- 
ergy shift from fullerenes to nanotubes can be explained 
by the influence of curvature on the coupling between a 
orbitals. In Fig. 2(b), we plot a detail of the low-energy 
region of the same spectrum, where the contribution of 
more neighboring cells has been plotted separately. The 
7r-plasmon peak is basically determined by coupling of 
the central cell with first neighbors, while the contribu- 
tion of more distant cells overcomes that of first neigh- 
bors below this peak. This observation suggests that the 
position of the 7r-plasmon feature depends very little on 
the long range structure of the sample, while the sub 
7r-plasmon features are more sensitive to geometry on a 
larger scale. 

Reed and Sarikaya have measured loss spectra corre- 
sponding to a single nanotube of about 1.2 nm in diame- 
ter using a 100-keV beam focused over a region of 0.2 nm 
in diameter. Their sample is similar to the one used in 
our theoretical simulations. The spectrum corresponding 
to the probe passing through the tube presents a broad 
resonance around 17 eV, while a faint signature of a low 
energy resonance is hidden by the elastic peak. For large 
impact parameters (several nanometers away from the 
tube axis), the low-energy loss peak becomes sharper and 
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the position of the a + tt plasmon feature is red shifted 
to 15 cV. Although the comparison of these data with 
those derived from our model is not straightforward be- 
cause of the different type of beam, one can claim good 
qualitative agreement between them. The contribution 
of electrons with large impact parameter can be selected 
by collecting the electrons scattered within a small scat- 
tering angle. In Fig. 3(a) we show three spectra corre- 
sponding to different values of the collector aperture 9 C . 
As the angle 9 C decreases, the ir peak becomes relatively 
higher, while the center of the a + it peak shifts down 
in energy. This dependence of the intensity of the dif- 
ferent peaks of the spectrum on the scattering angle has 
been also experimentally confirmed for Ceo- 5 A simple 
explanation based on the classical EELS theory has been 
proposed by Reed and Sarikaya, 14 namely, that the im- 
pact parameter dependence of the energy loss probability 
goes like 

P(o))-e~ 2 ^, d>v/uj, (25) 

where d is the distance relative to a given scattering 
center. 38 This means that distant collisions (i.e., small 
momentum transfers) excite low-energy modes more ef- 
ficiently, a fact that explains the positive dispersion of 
the loss peaks. The analysis of the 7r-plasmon feature of 
Fig. 3(b) also shows a weak positive dispersion, which is 
in qualitative agreement with momentum resolved EELS 
measurements in multishell fullerenes and peapods. 9 ' 15 

The EELS simulation corresponding to spherical 
fullerenes is more straightforward, since the system has 
now a finite number of C atoms. 39 In Fig. 4, we show the 
resulting computed loss spectrum in Cgo as compared to 
experimental data reported by Keller and Coplan, 5 mea- 
sured in the gas phase, so that, as discussed in Sec. I, our 
simulation based upon a broad beam configuration suits 
best the experiment. The energy of the electrons is 1 keV 
(v = 8.6 a.u) and the collector aperture is 1.5°, which 
is small enough to justify the neglect of recoil, even at 
such a relatively low velocities. Although the theoretical 
spectrum seems to be shifted up in energy, the simula- 
tion reproduces reasonably well the relative position of 
the peaks (including the small sub ir plasmon structure) 
as well as their relative intensity. In Fig. 5, we compare 
the results of our theory to the experiments of Kuzuo 
and coworkers 11 for crystalline Cgo- The probe energy is 
now 60 keV (v — 61 a.u.). We obtain good agreement in 
this case as well, suggesting that the coupling between 
neighboring target fullerenes is weak. This spectrum is 
very similar to those corresponding to C70, C76, and Cs4, 
with only small differences in the low energy region be- 
low the rr-plasmon feature, as shown in Fig. 6. The n 
plasmon around 6.7 eV is very similar in all cases, ex- 
cept for a small red shift as the number of atoms in the 
molecule increases, in agreement with experiments us- 
ing crystalline fullerenes, 11 where the ir peak lies in the 
6.4 — 6.1 eV interval. (Nevertheless, one has to keep in 
mind that the interaction between neighboring fullerenes 



in the solid phase can play a role that is not described 
by the present calculations.) The almost identical shape 
of the 7r-plasmon peak in this series (Fig. 6) is in con- 
trast with the strong differences at lower energy losses, 
as a consequence of the fact that the local atomic struc- 
ture of all of these molecules is very similar up to first 
neighbors (which give the dominant contribution to the ir 
plasmon), while differences in geometrical ordering over 
larger distances play a significant role at lower energies. 

In summary, we have obtained the screened interaction 
of a system of polarizable atoms to an external electric 
field and have applied this function to study the energy 
loss experienced by an electron (within a fully quantum- 
mechanical treatment of the probe) passing near carbon 
fullerenes and nanotubes. The spatial range of the cou- 
pling between dipoles explains the sensitivity of EELS 
to the structural details of the carbon nanostructures. 
The good agreement between theory and experiments 
is encouraging regarding the applicability of the derived 
screened interaction to other problems involving geomet- 
rically well-defined atomic structures. 
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FIG. 1. Polarizability of a carbon atom in graphite for the 
electric field perpendicular (a) and parallel (b) to the atomic 
planes. 
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FIG. 2. (a) Energy loss probability per carbon atom in a 
(10, 10) carbon nanotube. The primary electron energy is 100 
keV and the collector aperture half-angle is 9 C — 1 mrad. The 
solid curve represents the total energy loss probability, while 
the rest of the curves stand for the partial contribution of the 
central cell, the first neighboring cells, and the central cell 
plus the first neighboring cells (see labels and discussion in 
the text), (b) Detail of the low-energy region of the spec- 
trum in (a), showing separate contributions of up to the fifth 
neighboring cells. 
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FIG. 3. (a) Energy loss probability per carbon atom in 
a (10, 10) carbon nanotube for three different collector aper- 
ture half-angles (see labels) . The primary electron energy is 
100 keV. (b) Detail of the 7r-plasmon peak in (a) showing a 
positive dispersion with increasing aperture size. 



FIG. 5. Comparison between theory (solid curve) and ex- 
periment (circles, taken from Ref. [11]) for the energy loss 
probability of 60-keV electrons in solid Ceo- The collector 
aperture half-angle (not specified in the experimental refer- 
ence) has been set to 20°. 




FIG. 4. Comparison between theory (solid curve) and ex- 
periment (circles, taken from Ref. [5]) for the energy loss 
probability of 1-keV electrons passing through a Ceo sample 
in the gas. The collector aperture half-angle is 1.5°. 



FIG. 6. Electron energy loss spectra calculated for differ- 
ent spherical carbon molecules. The plots have been shifted 
upwards to allow comparison. The primary electron energy is 
60 keV. Electrons are collected up to a maximum scattering 
angle of 1 mrad. 
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